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A multi-body flight simulation for the Phoenix Mars Lander has been developed that 
includes high fidelity six degree-of-freedom rigid-body models for the parachute and lander 
system. The simulation provides attitude and rate history predictions of all bodies 
throughout the flight, as well as loads on each of the connecting lines. In so doing, a realistic 
behavior of the descending parachute/lander system dynamics can be simulated that allows 
assessment of the Phoenix descent performance and identification of potential sensitivities 
for landing. This simulation provides a complete end-to-end capability of modeling the entire 
entry, descent, and landing sequence for the mission. Time histories of the parachute and 
lander aerodynamic angles are presented. The response of the lander system to various wind 
models and wind shears is shown to be acceptable. Monte Carlo simulation results are also 
presented. 
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I. Introduction 

A multi-body flight simulation for the Phoenix Mars Lander has been developed that includes a high fidelity 
model for the parachute. During final descent, the Phoenix Lander relied heavily on its onboard radar to measure 
distance to the ground and vertical and horizontal velocities. The radar was sensitive to rotation rates of the antennae 
and was susceptible to a reduction in performance for angular rates greater than 100 deg/s. Use of the radar system 
places additional emphasis on accurate modeling of the parachute dynamics because the system depends critically 
on the exact orientation of the lander relative to the surface. To accurately model the dynamics of this system 
requires a high fidelity end-to-end simulation of the entire entry, descent and landing (EDL) sequence. Flight 
simulation of a vehicle with parachutes is mathematically complicated and not easily characterized using point-mass 
or single rigid body approaches. It involves multiple bodies, some of which are very flexible, flying in close 
proximity with significant interaction effects. Because the bodies fly so close to one another, the flowfield around 
each vehicle has an aerodynamic effect on the others. There are also direct, physical connections via lines. These 
lines exert a tension force that is a function of both the relative positions and the relative velocities of the connected 
bodies. 

This multi-body simulation of flight under a parachute was incorporated into a larger simulation of the entire 
entry, descent and landing (EDL) sequence. In the current work, the parachute, from the riser confluence point back, 
is treated as a rigid body, with the interaction forces between the parachute and other rigid bodies included. A 
similar treatment was used to simulate Viking drop-test flights 1 , as well as the Mars Exploration Rover (MER) entry 
trajectory. Detailed modeling of the dynamic aero- structural interaction of a flexible body such as a parachute is 
impractical in a full entry simulation. By treating the parachute as a rigid body, it is possible to create an end-to-end 
simulation from the last trajectory correction maneuver before atmospheric entry to touchdown. This simulation 
provides a realistic model of the descending parachute/lander system dynamics. In particular, position, attitudes and 
attitude rates are simulated and used to drive high fidelity models of the onboard radar and flight software. This 
exercise of the onboard system provides a method to assess system performance and to identify potential 
sensitivities for landing. Other issues such as recontact with jettisoned elements (heatshield, backshell, parachute 
mortar covers, etc.), design of parachute and attachment points, and desirable line properties can also be assessed 
using this simulation. 


II. Entry, Descent and Landing Sequence 

The Phoenix mission performed a direct entry at Mars. The vehicle entered the Martian atmosphere then 
decreased velocity via atmospheric drag until a Mars heritage disk-gap-band parachute was deployed. The entry, 
descent and landing sequence of events is shown in Fig. 1. 
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Figure 1. Phoenix EDL sequence of events. 
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The reference mission begins when the onboard entry state is initialized at E-10 min prior to atmospheric entry. 
Atmospheric entry is assumed to occur at a radius of 3522.2 km (radial altitude of 125 km). The entry inertial 
velocity is 5.6 km/s and the inertial flight path angle is -13.0 deg. When the vehicle is slowed to Mach 1.65 
(approximately 370 m/s), the parachute is ejected from the back of the spacecraft by a mortar. The mortar is sized to 
have an initial separation velocity relative to the spacecraft of 37.92 m/s. The parachute is packed in a bag (assumed 
massless for modeling analysis) until the lines connecting the parachute to the vehicle become taut, at which point, 
the bag is discarded and the parachute begins to inflate. 

Following deployment of the parachute, the heatshield is jettisioned and the lander legs are deployed. Once the 
lander legs are deployed, the radar is activated. The parachute/lander combination descends until lander separation, 
which occurs at a nominal altitude of 940 m. This separation altitude is determined by the onboard flight software, 
after which the lander performs a tip-up maneuver and then descends to the surface under power of its main descent 
engines. 

The Phoenix lander was actively controlled but unguided during its entry, descent and landing. The control 
system was tuned to minimize thruster firings prior to the deployment of the parachute. All of the critical events of 
the EDL phases were triggered by onboard flight controls. Reference 2 provides an overview on the development 
and evolution of the EDL sequence. 

This multi-body simulation modeled all the EDL phases of Phoenix: entry preparation, hypersonic, parachute, 
and terminal descent, including the onboard sequencing. Reference 3 provides a detailed description of many of the 
models incorporated into this simulation. However, the simulation described in Ref. 3 modeled the capsule and 
parachute as a single, rigid body, while the simulation described in this paper models the capsule and parachute as 
separate bodies (i.e., a multi-body simulation). Consequently, actual two-body capsule/parachute dynamic behavior 
is obtained allowing a more realistic determination of the performance characteristics of the radar for landing. While 
both the single body and multi-body simulations were used throughout design and operations of the mission, the 
single body was primarily used during the final days leading to Mars entry due to its faster run time. Performance of 
the radar was of great concern, since successful landing was so reliant on successful operation of the radar. The 
geometry of the lander with a supersonic parachute deployed is shown in Fig. 2. 
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Figure 2. Phoenix parachute configuration. 

III. Parachute Description and Modeling 

The Phoenix parachute is a disk-gap-band parachute based on the Viking design. It is 11.73 m in nominal 
diameter. Note that the projected diameter of the parachute in flight (as shown in Fig. 2) differs from the nominal 
parachute diameter. The aerodynamics of the parachute are represented as coefficients of normal and axial force and 
pitching moment about the confluence point of the parachute suspension lines. The normal and axial force 
coefficients and the pitching moment coefficient vary with total angle of attack. The normal and axial force 
coefficients of the parachute are shown in Figs. 3 and 4 as a function of total angle of attack, and the pitching 
moment is shown in Fig. 5. In figures 3, 4 and 5, upper and lower bounds are shown. These bounds represent the 
upper and lower bounds used in the Monte Carlo analysis discussed below. The dispersions for these three quantities 
were correlated so that when the normal force was perturbed in a positive direction the axial force and pitching 
moment were perturbed in a negative direction and vice versa. 
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Figure 3. Phoenix parachute normal force coefficient. 



Figure 4. Phoenix parachute axial force coefficient. 
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Figure 5. Phoenix parachute pitching moment coefficient. 



Time from Mortar Fire, sec 
Figure 6. Phoenix parachute inflation profile. 


Parachute inflation takes place over 1.2 seconds based on an inflation profile from Ref. 5. The inflation profile is 
shown in Fig. 6. This inflation ratio is used as a multiplier on each of the aerodynamic coefficients for the parachute, 
from the time of mortar fire until full inflation. The axial force coefficient was modified according to the Mach (M) 
variation shown in Fig 7. 



M 

Figure 7. Phoenix parachute Drag variation with Mach Number. 
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IV. Approach 

The underlying simulation software used for Phoenix is the Program to Optimize Simulated Trajectories II 
(POST II). POST II is the latest major upgrade to POST. 6 POST was originally developed for the Space Shuttle 
program to optimize ascent and reentry trajectories. Over the years, it has been upgraded and improved to include 
many new capabilities. POST II relies on most of the technical elements established by POST, but the executive 
structure has been reworked to take advantage of today’s faster computational capabilities. The new executive 
routines allow POST II to simulate multiple vehicles simultaneously, and to allow a combination of three degree of 
freedom (3DOF) vehicles with six degree of freedom (6DOF) vehicles within a single simulation. The multi-vehicle 
POST II capability was enhanced to allow interaction of forces and moments between the separate vehicles. As a 
result, parachutes interacting with other vehicles can be modeled allowing for a multi-body simulation capability. 
The lines connecting a vehicle and a parachute are modeled as massless spring-dampers. The springs can be attached 
at any point on the vehicle. No moments are applied except those due to force application away from the center of 
mass. Each line connects an attachment point on one vehicle to an attachment point on another vehicle and provides 
a tension-only force. When the lines are stretched, tension in the lines is determined as a function of strain and strain 
rate. Each of the spring-damper lines has an unstretched length and if the separation distance between the two 
attachment points is less than the unstretched length, the line tension is zero. Compressive damping forces are 
limited to be less than spring forces, so no line ever exerts a net compressive force. 

For the results shown in this report, damping and stiffness were both assumed to be linear. For the single riser, 
both subsonic and supersonic, the stiffness used was 60,000 N/m, while for the triple bridle legs (each) the stiffness 
was 47,000 N/m. The damping in the riser was 600 N/m/sec, while the damping in the triple bridle legs (each) was 
470 N/m/sec. The unstretched length of the riser was 2.778 m, while that of the triple bridle legs (each) was 0.914 m 
(see Fig. 2). The triple risers were attached to a plate located 1.7305 m aft of the capsule nose. Each riser was 
located 0.1 19 m from the axis of symmetry of the capsule and the three bridle legs were located symmetrically, each 
clocked 120 deg from each other. The canopy and suspension lines were modeled as a single rigid body in order to 
alleviate the extreme numerical sensitivity, which was evident if the lines were modeled individually. 

V. Model Verification 

In order to verify the multi-body model as implemented in POST II, a series of increasingly complex test cases 
was performed. These test cases were intended to prove that the POST II model is implemented correctly by 
evaluating its performance on problems that could be verified by other means. The test cases started with a simple 
vertical drop from rest of the fully deployed parachute and capsule, and then gradually increased the simulation 
complexity to include parachute opening, lander deployment, retro-rocket firing, and other effects. The POST II 
model had been previously compared and validated to MATLAB-based and another multi-degree-of-freedom 
simulations. 8,9 In each case the agreement between the simulations was excellent. The verification work is reported 
in more detail in Refs. 8 and 9. 


VI. Results 

After completion of the verification phase, the specific parachute configuration of the Phoenix Lander was 
incorporated into the simulation. The simulation was then used for a variety of analyses to characterize the overall 
performance of the Phoenix end-to-end entry, descent, and landing system. The response of Phoenix to a series of 
wind shears was investigated. The response at parachute deployment was simulated and compared to the response 
obtained using a simple “Wrist Mode” model of the parachute. The Wrist Mode model is a single-body model which 
treats the parachute as a pure drag device acting on the entry capsule. In addition, sensitivities to three different wind 
profile models were analyzed and the responses were assessed. Finally, the model was incorporated into a Monte 
Carlo simulation and the response was examined in a statistical sense. 

A. Wind Shear Response 

The behavior of the lander and parachute as they descend through a wind shear was examined. Figure 8 shows 
the wind profiles used in this test suite. For this test suite, the parachute and capsule were assumed to be initially 
falling straight down with a velocity of 60.6 m/s. Four cases were examined with wind shears beginning at 4, 3, 2 
and 1 km above ground level. The vertical angle response of the capsule is shown in Fig. 9 and the lateral rotation 
rate is shown in Fig. 10. It is clear from the angle and rate responses that the capsule is well damped for this type of 
excitation. In Fig. 9, it can be seen that after an initial transient attitude excursion of approximately 30 deg, the 
capsule settles to less than 10 deg within about 250 m. Similarly, Fig. 10 shows that transient rates of approximately 
40 deg/sec damp to less than 10 deg/sec within 250 m. 
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Figure 8. Wind shear profiles. 
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Figure 9. Capsule angular response to wind shear. 
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Figure 10. Capsule lateral rotation rate response to wind shear. 

B. Parachute Deployment and Inflation 

Since the opening of the parachute is a critical event, it was of particular interest. A set of simulations was 
performed to determine the vehicle response if the parachute took longer than expected to deploy and inflate or if the 
parachute "breathed" as it initially inflated. This "breathing" was modeled as a sinusoidal variation of the parachute 
drag coefficient about the nominal coefficient. These simulations included the nominal PEDLAM wind profile 
discussed below. The delayed opening case was modeled by stretching the inflation profile from Fig. 5 by a factor of 
five. Instead of opening in just over a second, the delayed opening took nearly 6 seconds. Figure 1 1 shows the drag 
force on the parachute for the cases considered. The delayed opening case lowered the opening loads slightly while 
the breathing cases raised the maximum drag force from approximately 47 kN to 48 kN. As a result, the capsule 
body rates are excited by the opening event, but they damp out very quickly as seen in Fig. 12. 
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Figure 12. Capsule rates after parachute deploy. 


Figure 13 shows the capsule total angle of attack from the time of parachute deploy to the ground. None of the 
off-nominal parachute openings had a significant effect on the angle of attack history. The attitude remains small 
(less than 30 deg) until lander separation at approximately 420 sec when the tip-up maneuver is initiated increasing 
the total angle of attack to -65 deg prior to the start of the gravity turn for landing. The altitude- velocity histories of 
these five cases are shown in Fig. 14. Immediately following parachute deploy, the delayed opening case, in 
particular, diverges from the nominal, but by about 9 km above ground level (AGL), all of the trajectories are 
essentially nominal, and all cases performed successful landings. 
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Figure 13. Capsule total angle of attack after 
parachute deploy. 



Relative Velocity, m/s 

Figure 14. Altitude-velocity histories for off- 
nominal parachute deployment. 


C. Wind Models 

Because the winds have such an impact on the behavior of the system while on the parachute, several wind 
models were investigated beyond the artificial wind shears described above. For Phoenix, three wind models were 
developed and examined: the Large Eddy Simulation (LES) wind model, the Turbulent Kinetic Energy (TKE) wind 
model and the Phoenix EDL Atmospheric Model (PEDLAM) wind model. The nominal north and east wind profiles 
for each of these models is shown in Fig. 15. It is apparent from Fig. 15 that the PEDLAM wind model has more 
high frequency content than the other two profiles. This model was the most stressful wind profile considered and 
was adopted as the reference for all subsequent analyses for conservatism. 

The total angle of attack for the capsule is shown versus time in Fig. 16. It is clear from Fig. 16 that the 
PEDLAM wind model has the most impact on the attitude of the capsule (again the tip-up maneuver is observed at 
-420 sec). However, the LES winds result in a largest angle of attack at touchdown. Note that the parachute trims 
near 10 deg total angle of attack. Since the aerodynamics of the parachute dominates the capsule/parachute 
combination, the parachute attitude is reflected in the trim attitude of the capsule as well. 
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Figure 15. Nominal wind model profiles. Figure 16. Capsule total angle of attack for various 

wind profiles 

Figures 17 and 18 show the pitch angle responses of the capsule and parachute respectively for each of the wind 
profiles. Comparison of Figs. 17 and Fig 18 shows that the parachute is much less susceptible to high frequency 
motion than the capsule beneath it, and that the PEDLAM wind excites both capsule and parachute more than the 
other winds. In each of the wind cases considered in this analysis, the vehicle landed well within requirements. 
Table 1 summarizes some of the relevant touchdown statistics from these wind profile sensitivities. The vertical 
touchdown velocity requirement is greater than 1.4 m/s and less than 3 m/s. The horizontal touchdown velocity limit 
is less than 1 m/s. It is clear that the vehicle makes a safe landing for each of the wind options. The TKE winds 
excited the largest attitude dispersions on the parachute while the LES winds excited the largest rates. 
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Figure 17. Capsule pitch angle wind comparison. 



Figure 18. Parachute pitch angle wind comparison. 


Table 1. Wind Comparison Statistics 



Wind 

Quantity 

None 

TKE 

LES 

PEDLAM 

Final Velocity (m/s) 

2.611 

2.837 

2.515 

2.802 

Final Horizontal Velocity (m/s) 

0.584 

0.504 

0.795 

0.369 

Final Vertical Velocity (m/s) 

2.545 

2.792 

2.385 

2.778 

Max Attitude from Nadir on Parachute (deg) 

17.072 

28.890 

17.222 

22.751 

Max Rate on Parachute (deg/s) 

77.461 

77.855 

93.703 

50.713 

Max Rate after Heatshield Separation (deg/s) 

20.481 

24.703 

30.833 

27.595 

Rate at Lander separation (deg/s) 

7.129 

4.682 

3.260 

6.071 
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D. Monte Carlo Simulation 

The parachute model was incorporated into a Monte Carlo simulation. The Monte Carlo simulation was very 
similar to the single-body analysis described in Ref. 3. The input dispersions were the same as those in Ref. 3, 
except for the inclusion of the multi-body dispersions (e.g., parachute aerodynamic data, line properties, etc.) 
described previously in Sections II through IV, along with the specific inputs shown in Table 2. 


Table 2. Monte Carlo Inputs 


Variable 

Distribution 

Units 

Mean 

3 Sigma/Limits 

Line Stiffness Multiplier 

Normal 

- 

1.0 

0.1 

Line Damping Multiplier 

Normal 

- 

1.0 

0.1 

Entrained Air Volume 

Uniform 

m 3 

602.203 

±403.145 

Parachute Ixx, Iyy, Izz multiplier 

Normal 

- 

1.0 

0.3 

Parachute Ixy, Iyz, Ixz values 

Normal 

Kg*m 2 

1.0 

30 

Parachute Inflation Rate multiplier 

Normal 

- 

1.0 

0.1 

Mortar Velocity 

Normal 

m/s 

35.65 

6.75 

Mortar Pitch Relative to Capsule 

Normal 

deg 

0.0 

2.5 

Mortar Pitch to Parachute 

Normal 

deg 

180 

2.5 

Mortar Yaw 

Normal 

deg 

0 

2.5 


Some selected relevant statistics from the multi-body Monte Carlo simulation are shown in Table 3. It is 
apparent from Table 3 that relatively small angles of attack at parachute deploy (<10 deg) led to fairly high rates 
through heatshield separation (> 100 deg/s). However, these high rates damped out very quickly. By the time of 
heatshield separation, the 99% high rates were less than 50 deg/s and by lander separation, the 99% high rates had 
reduced to less than 35 deg/s. These results showed that the multi-body parachute system reliably damps out 
rotational motion of the capsule. This rotational damping was of intense interest to the project because of the 
sensitivity of the radar to angular rates greater than 100 deg/s. This Monte Carlo resulted in 0.3 percent of cases 
violating constraints. The reasons for these violations were: horizontal velocity at touchdown too high (3 cases out 
of 2000), angle of attack too large at parachute deploy (2 cases out of 2000), and dynamic pressure too high at 
parachute deploy (1 case out of 2000). This was considered an acceptable risk to the project because of conservatism 
built into the simulation. This multi-body simulation provided sufficient confidence in the descent parachute/lander 
system dynamics to alleviate concerns of radar sensitivity or performance. The Phoenix Lander successfully landed 
on Mars on May 25, 2008. Details of the Phoenix EDL reconstruction and actual performance can be found in Ref. 
10 . 


Table 3. Monte Carlo Output Statistics 


Parameter 

Units 

99% Low 

Mean 

99% High 

Max Total AoA from Chute Deploy to Heatshield sep 

deg 

8.797 

15.870 

26.195 

Max Pitch Rate from Chute Deploy to Heatshield sep 

deg/s 

8.169 

44.944 

121.621 

Max Yaw Rate from Chute Deploy to Heatshield sep 

deg/s 

8.620 

38.513 

105.168 

Total AoA at Chute Deploy 

deg 

0.251 

2.348 

7.656 

Total AoA at Heatshield sep 

deg 

0.808 

7.881 

19.969 

Total AoA at Lander sep 

deg 

1.433 

11.681 

24.335 

Pitch Angle at Lander sep 

deg 

-88.367 

-77.58 -62.69 

RSS Pitch Yaw rate at Chute Deploy 

deg/s 

0.337 

4.349 

15.165 

RSS Pitch Yaw rate at Heatshield sep 

deg/s 

1.463 

15.267 

47.916 

RSS Pitch Yaw rate at Lander sep 

deg/s 

1.426 

11.741 

33.038 


VII. Conclusion 


A multiple rigid-body parachute simulation model for the Phoenix entry, descent, and landing has been 
developed. The model includes dynamics caused by interacting lines and has been incorporated into an end-to-end 
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simulation of the entry, descent and landing. This approach produced a higher fidelity simulation, especially with 
regard to capsule attitudes and rates while on the parachute. The capsule attitude and attitude rate were of particular 
interest to the mission team because of the sensitivity of the radar to large off-nadir angles and to large angular rates. 
It was seen that Phoenix was robust to wind shears of 20 m/s while descending on the parachute. It was also seen 
that several off-nominal parachute openings could be accommodated by the system. In particular, when the 
parachute opening was drawn out by a factor of five and when an oscillatory opening (breathing) was enforced, the 
system continued on to a safe landing. Several wind models were investigated. The PEDLAM wind model was the 
most stressful to the system, but it was not so stressful as to prevent safe landing. Finally, the model was used in a 
Monte Carlo simulation which showed that even very high rates soon after parachute deployment were quickly 
damped down to acceptable levels. 
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